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Abstract 



Q I The scaling form of the whole distribution P{D) of the random diffusion coefficient D{x) in a model 



of classically diffusing particles is investigated. The renormalization group approach above the lower 
critical dimension d = is applied to the distribution P{D) using the n-replica approach. In the 
annealed approximation (n = 1), the inverse gaussian distribution is found to be the stable one under 
rescaling. This identification is made based on symmetry arguments and subtle relations between this 
model and that of fiuctuating interfaces studied by Wallace and Zia. The renormalization-group fiow 
for the ratios between subsequent cumulants shows a regime of pure diffusion for small disorder, where 
P{D) — > 5{D — D), and a regime of strong disorder in which the cumulants grow infinitely large and the 
diffusion process is ill defined. The boundary between these two regimes is associated with an unstable 
fixed-point and a subdiffusive behavior < x^ >~ t^^'^''^. For the quenched (n — > 0) case we find that 



unphysical operators are generated raising doubts on the renormalizability of this model. ImpUcations 
to other random systems near their lower critical dimension are discussed. 



I. Introduction. 

The asymptotic behavior of classical particles diffusing in a disordered media has been the 
focus of many recent investigations [1]. The models which have attracted most attention are 
formulated in terms of random transition probabilities or local random quenched forces. In most 
of these problems rf = 2 is a special critical dimension and most attention has been foucused 
on the anomalous diffusion which occurs usually for d < 2 [2-7]. The anomalous diffusion is 
manifested by the long-time behavior of the form: 

{x\t))^e^, (1) 

with 2/x < l(2yU > 1), the behavior is said to be sub- (super-) diffusive. Attention has been 
devoted also to a simpler model in which the diffusion constant D{x) is a local quenched random 
variable and the diffusion equation for the local density p(x, t) is: 

^-P^ = V.[D{x)Vp{x,t)]. (2) 

This equation describes as well the continuum limit of a random resistor network in which 
D[x) is the local conductance. In this context one can define an Hamiltionian £^ = | / d'''xD[x){yV{x)Y, 
where V{x) is the local voltage. The physical problem in which the D[x) are themselves fluctu- 
ating dynamical variables required averaging the partition function (annealed average), while if 
the D{x) are frozen, it is the free-engery which should be averaged (quenched average). 

The local D{x) are expressed as a sum of a uniform part D and a deviations part AD{x): 

D{x) = D + AD{x) = D + 25D{x), (3) 

{SDix)SD{y)) = {6D')S{x-y), (4) 

(where 5D{x) = — j^ was introducted for computational convenience). This model was con- 
sidered somewhat trivial is because for weak disorder, -—- << 1, there is an asymptotic normal 



diffusive behavior at any dimension d > (the irrevlevance of the disorder near the pure diffusive 
fixed-point is exphcitly shown below). 

However, this does not preclude a non-trivial behavior for the strong disorder regime. The 
existence of such a regime will be manifested by a finite basin of attraction of the pure diffusive 
fixed-point in the space of all possible probability distribution P{D) for the locally uncorrelated 
diffusion constants. This space may be represented by the couplings g'^^'> related to the cumulants 
of the distribution by: 

{5D{x)5D{y))^ = 2\g^^^5{x = y), 
{SD{x)SD{y)W{z))^ = 3\{-l)g'^'^6{x = y = z), 



(5) 



{5Dix,)---6Dixk))^ = k\i-l)''g^'^5ix,=X2 = --- = x,). 



The central question we set ourselves to address in this work is whether there exists a non- 
trivial distribution that will be associated with the "critical" behavior on the separatrix board- 
ering this basin of attraction in the cumulants infinite-dimensional space. Of course there may 
be more than one such distribution or even a whole family of them. That depends on the 
renormalizat ion-group fiow on the separatrix ("critical manifold") itself. 

However, even the existence of such one scale invariant distribution is far from trivial since 
the theory may be either non-renormalizable (if an infinite number of new relevant couplings 
are generated) or may be renormalizable in a larger "unphysical" space (if relevant operators, 
distinct from the above cumulants, are generated). 

Similiar questions have arisen recently in the study of quantum diffusion (localization) where 
it was shown that the theory is renormalizable if the whole distribution of the conductance is 



considered [8-10]. The distribution takes three different forms depending on the system being 
metalhc, insulating, or at criticahty in between [8]. Similar questions of renormalizability always 
arise when a random system is considered near its lower critical dimension. For non-random 
models symmetry ensures the full renormalizibity near the lower critical dimension in which they 
have an infinite number of marginal operators mixed by the RG flow s, (e.g. 0{M) symmetry 
[11] for the M-component Heisenberg model in d = 2 + e). For random systems the effective 
replicated Lagrangian has a larger number of couplings without having larger symmetries to 
impose constraints on the renormalized couplings to be related to each other (as to preserve 
these symmetries). So the question of whether random systems are renormalizable near their 
lower critical dimension is much more delicate. While the quantum diffusion in rf = 2 + e is an 
example of a renormalizable theory [8-10] , the random-field 0{M) model in d = A + e provides 
a counterexample [12]. 

We therefore hope that our systematic study of another such system will help to shed more 
light on this puzzling question. 

The paper is organized as follows: In the next chapter (II) the n-replicated Lagrangian will 
be derived. From this Lagrangian, which has the cumulants of the initial distribution related 
to its bare couplings, all correlation functions may be derived. In chapter III we analyze the 
so-called annealed approximation (n = 1) first by the standard RG approach and then by a 
search of the universal distribution based on symmetry. Chapter IV is devoted to the study 
of the quenched average {n -^ 0). The problems encountered are discussed. In chapter V we 
summarize the implications of these investigations for the diffusion problem in particular and for 
random systems in general. 



II. The replicated Lagrangian. 

Without losing generality, we assume all the particles to be at i^ = at the initial time t = 0, 
namely: 

p{x,0)=6{x). (6) 

The Laplace transform of p{x, t) is defined as: 



P{x,m') = / dtp{x,t)e'"' \ (7) 

Jo 

p{x, w?) obeys the equation: 

w?p[x, m^) — p{x, 0) = V • [D{x) vecV]p{x, m^), (8) 



or: 



my{x, m^) - V ■ [D{x)V]p{x, m^) = S{x). (9) 

We may identify p{x, m^) with the Green's function 

p{x,m^) = G{x,m^), (10) 

which is the Laplace transform of the time domain Green's function 

G(f,t) = (p(f,t)p(0,0)). (11) 

All scaling properties may be extracted from the large x and t (or small m?) behavior of these 
Green's functions. To be able to average these correlation functions over the realizations of the 
disorder we first express them as a functional integral in the standard form: 



G(f m^) = I'^MmO)e-'^'^^ .^2. 



where: 

d'x<P['^-V-{^V)]<P. (13) 

or, after integration by parts and neglect of an unimportant boundary term: 



m = jd'x(^<^'+^iym- (14) 

The averages over the disorder may be obtained by utihzing the rephca trick . The n- 
rephcated partition function is obtained by including n fields 0"(x), (a = 1,2, ....,n): 

Z"= /'n^0"(^)e"'^f'^"', (15) 
with: 

£"r] = /rfME[^(r)^ + ^(v0")^]}, (16) 



Since lim„_^o( -2^"") disorder = 1, we may eliminate the denominator in Eq.(|T2D, and calculate the 
average Green's function from: 

{G{x, m^))d^sorder = Hm / prod^P^ (f )0i(f )0^(O)e- ^'^'^"['^'^l (17) 

Since the quenched average is usually difficult to perform the "annealed" approximation 
is often used. It consists in averaging independently the numerator and the denominator in 
Eq.(^). This is also equivalent to keeping a single replica or averaging Z" in Eq. (1T5|) with 
n = 1. Therefore we shall, in the rest of this chapter, analyze the generalized field threories with 
£"[(/)"] for any n. The results of specific cases of annealed (n = 1) and quenched (n -^ 0) average 
are discussed in the next two chapters. Let us use the vector notation (p = (0^, 0^, ■ ■ ■■, 0"). The 
replicated Lagrangian may be separated into a free part and an interacting part, C^ = Co + Cm 
where: 

Co = lld''x{D{V$y + m\^r}, (18) 
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and: 

C,n = j d''x{5D{x){V^f}. (19) 

The average over the disorder yields: 

(e-^-)..We. = e/^^^i:r^.^'^'[(^^>^ (20) 

where g^^'^ are the couphngs defined in Eq.(^). We therefore have to apply the RG analysis to 
the following averaged partition function: 

{Z-)d^sorder = / X{Vr{S)e^v{- / rf'^x(- [m^^^ + D( V^)^] - Y. 9^'^"^ <^ff)} , (21) 

for simplicity we shall choose Z) = 1 in the forthcoming calculations. Also in the renormalization 
scheme we choose to keep this term constant. To make it dimensionless should carry dimension 
of [(/)] = L^-i . Hence [(V^)^] = L"'^, and [^(^)] = U^-^'^'^. 

Therefore under rescaling L — > Ljh , g^^^ -^ i,-ik-i)dg{k) ^^^^ g{k) ^^^ irrelevant near the free 
gaussian theory for any dimension d > 0. We also observe that d = is the critical dimension at 
which these couplings become marginal. The natural small parameter will thus be e = c?. In the 
next chapters we go beyond the dimensional analysis by utilizing first the RG approach and then 
an alternate approach based on possible symmetries of the non-trivial fixed-point distribution. 



III. The annealed (n=l) approximation. 

We begin our analysis with the consideration of the simpler annealed approximation. That 
corresponds to single component field 0. We begin by demonstrating the problems that arise 
in taking the RG route and then show how these difficulties may be circumvented to find the 
scale-invariant distribution. 
1. Renormalization-group approach. 



To go beyond the naive dimensional analysis we need to expand e^~ ^^^ ^''^'^' ^ -f in a power 
series, separate and integrate momenta A/6 < g < A (where A = 1/a is the boundary of the 
Brillouin zone) using e~^°, reexponentiate, and rescale all momenta q ^ hq as to obtain a new 
effective Lagrangian. This can be done diagrammatically and here we only give the recursion 
relations for /c = 2, 3 to order one loop (/ = Inh): 



dl 



= -eg^'> + 12g^') + 56ig^'>y, (22) 

^^ 2eg^^^ + 2A0g^'^^g^^^ + 24(7^ + AlQig^^^f, (23) 

Similar expressions may be derived for all g^''\ They will take the form: 

^ = Pi'\g^'\g^'\---,g^'^'^) = -{k-l)eg(>^^+(3l,'\g('\g^'\...,g('^^'^), (24) 

where Pq is the same beta function calculated at e = 0. Since Pi''' have both positive and 
negative terms each of them has zeros. The challenge is to find at least one common zero, 
gi'') = gi'')* for all k, where they all vanish simultaneously. We also note that the g^s will be 
generated under the RG recursions once the bare value of one of them is not zero. However, 
different couphngs (not cumulants of P{D) ) will not be generated and the space of all g^s is closed 
under the RG transformation. These equations also tell us that at the fixed-point g'^''^* ~ e'^^^, 
and we may express them in terms of coefficients at such that: 

^W* = ake'-\ (25) 

Using the equations /3^(^^ = all a^ may, iteratively , be related to 02, e.g. : 

03 = 7^(1-5602) 

^3 _, „ 52 o 

04 = 100302 a2, etc . 



We can also deduce some information on the behavior away from the fixed-point. Suppose 
we try to rescale all g^''^ multiplicatively: 

^W = C(/)'=-^^('=)*, (26) 

very close to the fixed-point ({0) ~ 1, then from Eq.(p6D : 



di -^ dl "^ ^^ ^'^^^^ dV 



(27) 



But on the other hand from Eq.( 

-e{k - 1)^W*C'-' - C'/?^!^*} = C'-'(C -l){k- l)eg^^>. (28) 



dg _n. ^\Jk)*^k-l_^kn{k), 



dl 

Equating the two equations (^) and (pHj) above we obtain: 

^ = (C - 1)<, (29) 

from which we identify C as a relevant scaling field near the fixed-point [C, = 1) with scaling 

exponent </)^ = e. 

2. Invariant cumulant generating function: A search by symmetry. 

Instead of looking for the invariant distribution based on the recursion relations (which seems 
pretty hopeless) we shall base our search on global symm etry considerations. We shall look first 
at the moment generating function: 

oo 1 ^r^^J ^ Ttk 



f{u) = \og jyDP{D)e~^- = E(-l)' ^""""^If ^' u\ (30) 

Comparing with Eq.(^, we have: 

/(O) = 0, (31) 

/'(O) = 1, (32) 
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"^'^^=0 = 9^'^^^ (33) 



If we can find a function f{u) such that for 2u = (V0)^ the in tegral 



mm]-f^unmh (34) 



will be invariant under renormalization, then the corresponding P{D) (to be obtained by inverse 
Laplace transform of e-'^^"-' ) will also be invariant . 

So the problem has been reduced to finding the function f{u) for which £(/) is invariant. 
This search must be based on symmetry: one has to look for a function f{u) for which the 
integral in Eq. (|3^) is an invariant un der a symmetry operation which itself is preserved under 
the RG iterations (this insight comes from what is known about the role of the symmetry for 
RG near the lower critical dimension [11]). Since is a scalar this cannot be strictly an internal 
symmetry. It should be a symmetry that mixes the order parameter 0(af) and the coordinates 
X. The simplest one is to add as the {d + l)th component of a new vector : 

a;^ = (f,0(f)). (35) 

The simplest invariant is then just the total arc length in this space: 

— * 

S[f{^^}] =Jd'x^l + {Vphty, (36) 

which is invariant under rotations in the d + 1 -space, and is therefore a natu ral candidate for 
£[/]inEq.(|4D. 

This type of actions has been introduced and studied by Wallace and Zia [13] to model 
interfacial fluctuation of Ising systems in ci + 1 dimension (in this picture 0(x) is the height of 
the interface, it carries dimensions of length, and the overall action is made dimensionless by 
giving dimension of L^ to the temperature). 

11 



The inverse Laplace transform of f{u) ~ e^"*^^"'"^"'' ' is the inverse gaussian distribution [14]. 
This distribution has two free parameters which may be related to its average and its variance, 
in terms of which the universal renormalized distribution is : 

The form of this distribution remains invariant under the RG flow. The only parameter which 
changes is the ratio: 

Higher cumulants are expressed in terms of R and D : 

{AD''), = (2n - 3)!!/?"-^^". (39) 

The average I) is a redundant parameter which can be given any value D > 0. 

On the renormalized trajectory the RG flow for R are determined by the equation: 

which has -R* = e as its fixed-point. Thus the fixed-point distribution is (choosing D = 1 ): 

P*(D) = {^^f^ exp{-- ^^~^^' }, (41) 

'^ ' ^2'KD^e ^^ 2t D ^' ^ ^ 

this distribution with e = 1 is plotted in Fig.l. 

For R < e, the RG flow will be to smaller R or < AD^ >^ while keeping D fixed. The 
flow takes th e distribution to: 

lim( ^, )^/^exp{-— fraciD - DfDD} = 6{D - D), (42) 

as expected. 
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For R> e , the RG flow will be towards larger R . The moments < AD'' > diverage as i?^~^ 
and the diffusion process is ill-defined. (As R ^ oo, P{D) ~ D"^/^ for R~^ « D « R). 

It is also straightforward to identify ({I) disscussed above as the scaling field ({I) = R{l)/R*. 
Hence the crossover exponent for AR ~ i? — i?* is e . This is the only relevant direction near 
the fixed-point [15] and any variation from the fixed-point along other directions is irrelevant, 
namely the fixed-point is stable in all other directions in this parameter space and the critical 
manifold (the separatrix) has codimension one. Operators which break the -^ —(f) symmetry or 
which contain power of 0, rather than of its gradient, are also relevant [15] but are not important 
in the present context (besides m^0^ which is discussed next). 
3. Anomalous diffusion on the critical manifold. 

For any distribution on the critical manifold the asymptotic scaling behavior will be deter- 
mined by the RG flows near the flxed-point. These are expected to be different than the simple 
diffusion which occurs near the free gaussian flxed-point. 

At the flxed-point the fleld 0^ will acquire anomalous dimension in order to keep D (and all 
the rest of the couplings) flxed. The anomalous dimension will be exactly that associated with 
scaling fleld C~^(0 ^^ ^^e flxed-point may be interpreted as the renormalization factor which 
multiplies 0^ to keep the couplings flxed. Since ({I) ~ e^' = 6^, the rescaling of 0^ will be : 

02 _^ ^-102 _ ^d-2-.02_ (43) 

That induces an anomalous rescaling of the Laplace transform parameter m^ : 

m^ -^ m%'^+% (44) 

which by its deflnition scales as 1/t. Hence the time t will rescale as: 

t -> tb^^'% (45) 

13 



while length rescales as 6 ^, and therefore: 

(x^) ~r2~t2T7, (46) 

Hence the diffusion is anomalous with exponent fi in Eq.(|I|) given by: 

Although locally stable we cannot prove that this fixed-point is unique. This question could 
be explored looking for all possible functions /((V0)^) which are scale invariant. One way to 
approach the problem is to ask what are the non-linear symmetries preserved under RG. The 
"Euclidean" symmetry is one of them and the action in Eq.(36) is its only invariant. We cannot 
rule out, however, the possibility of other preserved symmetries with other invariants. It will be 
very interesting if this general question could be investigated more systematically. If other fixed- 
points exist they may not follow the simple gap scaling above, and a "multifractal" behavior 
cannot be ruled out. 

IV. The quenched average (n ^ 0). 

Our goal is to pursue the same succesful route we followed to the solution in the annealed 
approximation, for the quenched case as well. We therefore repeat the RG calculation for general 
n (n 7^ 1) to be analytically continued to zero. Here, however, we encouter a new difficulty: The 
space spanned by the couplings g^'^^ is not closed under renormalization. 

We shall explain it for terms to order (50)^ but similar behavior takes place for higher powers 
as well (see Appendix). The term in the bare Lagrangian is of the form: 

n d n d 

o=l i=l a,l3=l i,j=l 

14 



To order {g^'^''^^Y the contraction of two such terms give rise to a renormahzation of g^'^''^^ itself 
but also to the generation of term of the form: 

g(0,2)J2J2(^d^^ag^^P^^d^^ag^^^^^ (48) 

a,l3 i,j 

The complete recursion relations are given in the Appendix. Clearly such a term takes us outside 
of the realm of the diffusion model we began with. If we insist on still following the RG flow 
in the a larger number of unphys ical parameter space we can generalize the approach based on 
symmetry discussed above for the case n = 1 to general n, by looking at a vector in the d + n 
dimension space {xi, ■ ■ ■ , Xd, 4>^{x), ■ ■ ■ , 0"'(x)). The invariant area in this space is: 

S[<P''{x)]= [ d^x{detgf'^, (49) 



with(7,, = (5,, + E^=i-9.0"9,0-. 

Although no distribution may correspond to this function (since it cannot be a cumulant 
generating function) there is a non-trivial fixed-point for finite n. The incomplete one-loop 
analysis of Lowe and Wallace [16] yields R* = e/n and will diverges in the n ^ limit. If 
this behavior will surive for a larger numbers of loops, it would imply that e ven in the larger 
space the trivial fixed-point is the only one accessible by the RG approach (which still leaves the 
possibilty of a non-perturbative strong coupling behavior) . 

V. Conclusions. 

In view of the new understanding acquired in the quantum diffusion problem on the im- 
portance of renormalizing the full conductance distribution [8-10] (rather than the first two 
moments alone), we have addressed the question of the distribution of the diffusion constant in 
random classical diffusion. We have chosen here the simplest model for diffusion in disordered 
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media which has rf = as its critical dimension. Our analysis was performed on the replicated 
Lagrangian in which the cumulants of the distribution are related to the coupling constants. 

We first looked at the annealed approximation for which we could identify the invariant dis- 
tribution. We have found a non-trivial fixed-point along the renormalized trajectory separating 
a free gaussian diffusion fixed-point for small '^'^ ^ from a regime where the cumulants ratio 
divgerges < AD^''^^^ > / < AD'^''^ >^ oo. At the fixed-point itself an anomalous diffusion 
< x^ >~ t^~2 was found. 

For the quenched average n — ;► our results are so far negative: the theory is not renormaliz- 
able because unphysical terms are generated under renormalization. Such a behavior may occur 
in other random systems near their lower critical dimension. That may also indicate that an 
even more general approach (for example including the possibilty of replica symmetry breaking) 
may be necessary. Other investigations of the distributions of the appropriate physical quantities 
in other random models may shed more light and will be most worthy. 
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Appendix 

In this appendix, we present the recursion relations of two lowest order couplings in Eq. (|i9| ) 
for general n approach. As mentioned in chapter IV, we need to consider other (cross) terms 
like the one in Eq.(^) which does not appear in the original Lagrangian in Eq.jpTI). Including 
the terms generated, one should rewrite the Lagrangian as follows: 

£[0-(f)] = j d'x]- ^ d,rd,r - E E[^^''°H5^0"5.r^./5./) + ^(^'^n^ir^^/^.r^./)] 

- E E[^^''°H^^<^"'9,r)(5,/5,0^)(part^a/fc(/)^afc</)^) 

Q:,/3,7 i,j,k 

+^(°'^)(9i0"(9i0^)(o»j0''(9j0")(9fc/9fc0^)] + ■ ■ ■ higher power terms. (50) 

By diagrammatic calculations, we found the following recursion relations for g^"^'^") and g^^'"^"* ( to 
order one loop): 

^^ = _e^(2,o) + ^2g^m + (2^ + 2)^(i'2) ^ 3^(0,3) ^ 4o(^(2,o))2 ^ ^^g^ ^ 32)(?(°'2)(^(2,o) 

+ (2n + 14)((7(°'2))2^ (5^) 

^^ = -e^(o.2) + 8^(1,2) + (3^ + 6)^(0,3) ^ (i2n + 28)((7(°''))' + 64(?(0'2)(^(2,o) _^ i6(^(2-0))2_ (52) 

For the case n = 1, the distinction between the two terms in the Eq.(??) with coefficients 
(^(2,0)^ g(o,2) (disappears, so one can identify gf*-^-* in Eq.(??) as the sum of (7'-^'°-' and (y(('^'2) in Eq. (pOD . 
For the s ame reason, g'^'^'> in Eq.(^) is the sum of g'^'^'^\ g^^'"^^ and g^^'^\ Therefore the sum of 
the two equations above is ide ntical to Eq.(l2^), for n = 1. 
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Figure Captions 

Fig.l -The form of the invariant inverse gaussian distribution (Eq.(p7D) wit h parameters D = 1 
and < AD^ >= 1 (it is the critical distribution for d = 1 within the an nealed approximation). 
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